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Abstract. We revisit the long time dynamics of the spherical fully connected p = 2- 
spin glass model when the number of spins N is large but finite. At T = 0 where the 
system is in a (trivial) spin-glass phase, and on long time scale t > we show 

that the behavior of physical observables, like the energy, correlation and response 
functions, is controlled by the density of near-extreme eigenvalues at the edge of the 
spectrum of the coupling matrix J, and are thus non self-averaging. We show that the 
late time decay of these observables, once averaged over the disorder, is controlled by 
new universal exponents which we compute exactly. 
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1 . Introduction and model 

Non-equilibrium dynamics of spin-glass models has attracted much interest, both 
theoretically and experimentally, during the last 40 years [1]. In particular, the low 
temperature relaxational dynamics of such systems following a quench from a high- 
temperature conhguration is not only extremely slow but also displays “aging effects”. 
This means that the response of the system (as well as temporal correlations) depend 
strongly on the history of the sample since the temperature quench. It is useful and 
now customary to characterize quantitatively such aging effects by studying two-time t, t' 
observables, including in particular the local response R{t, t') and the auto-correlation 
function C{t,t') (see below for a precise dehnition) [2]. 

The analytical computation of these observables for hnite-dimensional spin glasses, 
and more generally for disordered systems in hnite dimension, is a very hard task and 
consequently there exist very few exact results for such models (see, e.g.. Ref. [3] in 
the context of disordered elastic systems). However, it was demonstrated that it is 
already instructive to study fully-connected (mean-held) spin-glasses, whose dynamical 
responses and correlations were shown to exhibit quite rich structures, similar to some 
extent to the ones observed in structural glasses for instance [1]. The purpose of the 
present paper is to revisit the late time dynamics of the simplest model of that type, 
namely the spherical two-spin model of large but hnite number of spins, in light of recent 
results obtained in the literature on random matrix theory (RMT). 

The spherical two-spin model, which we focus on here, is dehned by the Hamiltonian 

~ JjjSjSj , (1) 

where Sj, with z = 1, • • • ,N, are continuous spin variables constrained such that 

N 

E^ p) 

i=l 

and where J is a symmetric random matrix belonging to the Gaussian Orthogonal 
Ensemble (GOE) of random matrices (corresponding to the Dyson index /3 = 1): 
its elements are independently distributed Gaussian random variables with zero mean 
and variance proportional to 1/N (with this choice, the model ([^ has a meaningful 
thermodynamic limit N —)■ cx)). The model was originally introduced in Ref. |1] and 
studied by many authors ever since, see e.g. chap. 4 of the book [5]. Although it 
was shown to exhibit a phase transition at the critical temperature = 1 from a 
paramagnetic phase into a low temperature “spin-glass” phase the latter turns out to 
be in fact a “ferromagnetic in disguise” [5]. Namely, the system possesses only two 
ground states related by the symmetry s* —)■ —s* and the calculation of the free energy 
for this model using replicas does not involve a replica symmetry breaking which is a 
hallmark of the true spin-glass thermodynamics. As the minimum of the quadratic form 
([^ on a sphere is obviously given by (minus one half of) the largest eigenvalue Amax 
of the matrix J, the statistics of the ground state are governed by the Tracy-Widom 
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distribution for GOE, which describes the fluctuations of the largest/smallest eigenvalue 
in that ensemble |7]. Besides, it was recently proved that for T < T^, the fluctuations 
of the free energy are also given by the Tracy-Widom distribution, see [6]. 

Throughout this paper, we will be interested in the limiting case T = 0, 
which already contains many interesting aspects. Although, as explained above, the 
thermodynamics of the model in the low-temperature phase is too simple for a bona 
hde spin-glass, the corresponding dynamics is rich and has features of aging PEiiiniE]. 
That richness is attributed to a relatively non-trivial energy landscape topology due to 
the presence of many saddle points in the landscape with different indices k (the number 
of unstable directions), each associated with the eigenvalues A* < A^ax of the matrix 
J. Note that the presence of a magnetic held in the system leads to a simplihcation 
of the associated energy landscape, which gradually washes out the complexity of the 
relaxational dynamics (see Refs. [121 IlHl EH EH] for this and related questions). Our 
studies of the dynamics of the model ([^ will be intimately related to the setting of 
the original work by Cugliandolo and Dean [8] referred to as CD in the following. In 
particular, we will use the same convention for the distribution of the matrix elements 
Jij such that when N ^ oo, the distribution of the eigenvalues of the random matrix J 
is given by the Wigner semi-circle law on the bounded interval [—2, 2]: 

p(X) = . (3) 

To study the relaxational dynamics of this model ([^, it is convenient to diagonalize 
the coupling matrix J and to work with the time dependent projections of the spin 
conhguration s{t) = {si(f)}i<i<Ar onto the eigenvectors of J, which are denoted by 
sa(^), a belonging to the spectrum of J, A G Sp(J). The dynamics of the model is 
then dehned via a Langevin equation, which when projected onto the eigenvectors of J, 
yields a Langevin equation for the projections s\{t) that reads [f] 

^ = (A - zit))s^(t) + ftVi) + m , (4) 

where h\{t) represents an (inhnitesimal) external magnetic held - which is used here to 
compute the response function - and z{t) is a Lagrange multiplier which enforces the 
spherical constraint (|^. In Eq. (|^, ^\{t) is a Gaussian white noise of zero mean and 
correlations 

{Ut)^At’)) = 2T6,,y6{t-L), (5) 

where T is the temperature and where (■ ■ ■) denotes an average over the thermal noise. 

2. Summary of main results of the paper 

Here we consider the zero temperature dynamics, at T = 0, where the system is quenched 
at initial time f = 0 from a high temperature conhguration, described by a uniform initial 

f note that Eq. Q also describes the dynamics of a spherical height model, the Arcetri model recently 
introduced in m, in one spatial dimension. 
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condition 


Sa( 0) = 1 , VA e Sp(J) , 


( 6 ) 


which was shown to lead to a rich dynamical behavior [S]. In the large time limit, 
it is easy to see that the Langevin dynamics Q at T = 0, i.e. ^A(i) = 0 (assnming 
h\(t) = 0), will drive the system to relax to the configuration with minimal energy per 
spin (energy density) e(t —)■ cx)) = — |Amax where 


Amax = max A , 

AeSp{J) 


( 7 ) 


and correspondingly s\(t —>■ oo) —)■ The question that we will address in 

the present paper is the following: how does the system approach this final state of 
minimal energy? For instance, how does the average energy density e{t) approach its 
limiting value? To characterize the non-equilibrium dynamics of the system (|^ it is 
also useful to consider quantities depending on two times t, t'. Here we will focus on the 
zero-temperature limit of the spin-spin correlation function C{t,t') and on the response 
function R{t,d), t > t'\ 


C{ht') 

R{t,t') 


1 

N 

1 

N 


i=l A 

dsi{t) ^ 1 ds\{t) 

^ 5hAt') ~N^ 6hx(t') 

i=l ^ h=0 A ^ h=( 


( 8 ) 


and their corresponding disorder averaged values C{t,t') and R{t,t') where denotes 
an average over the disorder realizations, i.e. over the GOE random coupling matrix J. 
In Eq. (|^, and in the following, 'J2x denotes a sum over all the eigenvalues belonging 
to the spectrum of J. 

Intuitively, one expects that at large time, the dynamics in Eq. (|^ will be 
dominated by the near-extreme eigenvalues, i.e. the eigenvalues of the coupling matrix J 
that are close to A ma x laini. More precisely, we will see that the observables mentioned 
above (energy density, response and correlation functions) can be written in terms of the 
density of eigenvalues “seen” from Amax, the so called density of states (DOS) pBosi^, N) 
defined as [THl [TU] 

PDOs(r, N) = - A - r) . (9) 

The DOS ([^ was recently studied in detail for matrices belonging to the Gaussian 
Unitary Ensemble (GUE, corresponding to the Dyson index f3 = 2) [19], using semi- 
classical orthogonal polynomials, as well as for more general Gaussian /3-ensembles [20] , 
using mainly scaling arguments. In particular, it was shown that the behavior of the 
average DOS Pdos(?", N) exhibits two distinct behaviors depending on whether r ~ 0(1), 
or r ~ In the former case Pdos(^, N) is simply given by a shifted Wigner semi¬ 

circle law (see also im ), whereas in the latter Pdos (^"5 N) takes a non-trivial scaling form 
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reflecting the fluctuations at the edge of the Wigner semi-circle. This can be summarized 
as follows [IHl [20] 

_ f Pbuik(r’) = , r ~ 0(1) , 

Pdos(’", ~ \ (10) 

[ iV-V3pedge(A^2/3^) ^ ^ . 

The scaling function Pedge(3^) is presently known exactly only for complex Hermitian 
random matrices belonging to GUE [19]. For general Gaussian /1-ensembles, including 
the case /I = 1 mostly relevant for the spin-glass model studied here ([^, the exact 
evaluation of Pedge(^) remains an open problem. At the same time the asymptotic 
behaviors of such a scaling function was derived in 1201 and is given by 


ag X' 


^ , X —)■ 0 


X) 


( 11 ) 


-\/X ,x^oo. 

TT V ‘ 

The small-x asymptotics in ( [IT| is controlled by the repulsion between two neighboring 
eigenvalues, and the exact value of the amplitude is known only for (d = 2 where 
02 = 1/2. On the other hand, the large argument behavior can be obtained by matching 
the edge regime with the bulk regime given by the shifted Wigner semi-circle. 

In Ref. [S] (see also Ref. mm for a rigorous proof of these results), the authors 
studied the dynamics of the spherical 2-spin model in the thermodynamic limit N ^ oo. 
Most importantly the limit of large times t was studied after performing the limit 
N ^ oo. In this order of limits the large time behavior of physical observables can 
be obtained by replacing Pdos(’") dV) (which is self-averaging for r ~ 0(1)) by a shifted 
Wigner semi-circle pbuik(’") given by the first line of (10). Indeed, in the limit N ^ oo 
the contribution from the edge regime turns out to be negligible. In particular, it was 
shown that the average energy density e{t) approaches its limiting value algebraically 
as 


+ 2-^” 


( 12 ) 


for large time t. Such algebraic decay reflects the complexity of the energy 
landscape, with infinitely many saddle-points being operative in trapping the descending 
trajectories for a long time. Such saddle points typically have a “mesoscopic” index, 
that is the number k of unstable directions satisfy 1 k N, corresponding to 
eigenvalues of the matrix J mesoscopically far from the spectral edge. We notice that 
a similar decay oc 1/t, as in Eq. ( jl^ , was also observed in numerical simulations of the 
Sherrington-Kirkpatrick model [2T] (at T = Tf). 

In this paper, we are interested in clarifying what happens for such a system of large 
but hnite size N. In this case, we actually demonstrate that there exists a crossover 
time scale fcross ~ which separates two distinct relaxation regimes: the GD 

regime 1 as described above for the system in the thermodynamic limit N ^ oo and 
taking place for times t tcross from the later time regime 2 operative for t 3> Gross- 






Large time zero temperature dynamics of the spherical 2-spin glass model of finite sized 





Figure 1. Plot of N^^^Ae{t) as a function of for different values of iV = 

50,100,200 and 400. The data have been obtained by evaluating numerically the exact 
formula in Eqs. (23) and (30), and the averages have been performed by sampling 200 
independent GOE random matrices. The collapse of the curves for different values 
of on a single master curve (for t sufficiently large) is in good agreement with our 
analytical predictions in Eqs. (13) and (14). 


This latter regime is controlled by the small argument behavior of Pedge(2:) in Eq- 0. 
yielding in particular new exponents which we compute here. At this later stage the 
relaxation is dominated by just a few saddle-points in the landscape close in energy to 
the ground state and having of order of one unstable directions. 

For t ~ single time quantities like the average energy e{t) are characterized 

by a scaling function of the scaling variable which interpolates between these two 

regimes. For instance, we show that average excess energy Ae(t) = e(t)-|-Amax/2 behaves 
as 


ei{t) , t < Ar2/3 ^ 


Ae(t) ~ < 




iV2/3 


, t > 


(13) 


The function ei{t) was computed by CD in Ref. [H] and has the asymptotic behavior, 
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Figure 2. The three different regimes of the temporal behavior of the response 
function in the plane (note that the region t < t' is forbidden as a 

consequence of causality). The regime I corresponds to the regime studied by CD in 
[5] while the two regimes had not been studied before. In particular, in the regime III 
where the dynamics is stationary, the response function is completely determined by 
the density of near extreme eigenvalues of the coupling matrix J [see Eq. (15) ]. The 
dot indicates the region where the response function takes the scaling form given in 
Eq. (19) which interpolates between the three regimes. 


for large t, given in Eq. (12) while S{x) has the following asymptotic behaviors 

3 


£{x) 




A 

— , X -)■ CX) 
X-^ 


(14) 


where A is an (unknown) numerical constant. The regime x —)■ 0 naturally matches 
the CD regime in Eq. (12), while the other regime for x —?• oo, characterized by a 
different exponent, describes the late time behavior for a system of large but finite size. 
In Fig. we show a plot of iV^/^Ae(f), evaluated numerically from the exact formula 
derived in Eqs. (23) and (30) below, as a function of t/N'^/^ for different values of A, 
which conhrms our analytical predictions in Eqs. (13) and (14). It is important to stress 
that while in the regime 1, for t the quantities we study are self-averaging, we 

will argue that it is not the case for the late regime 2, for t > 


Similarly, two-time observables like the response R{t, t') and the correlation C{t, t') 
display different asymptotic behaviors in the three regions I, II and III in the (t, t') plane 
depicted in Fig. The hrst region I for times 1 t,T studied by CD [8], the 

region II for 1 <^t' t and hnally the region III for t, t’. In particular, 

in the region III the response function is stationary, i.e. it only depends on the time 
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difference t — t\ and R{t, t') is given by 

poo 

R{t, t') ^ Riii{t -t‘)= / dr , 


(15) 


which shows the physical relevance of the DOS in this problem. The scaling form for 


the DOS in Eq. (10) induces the following scaling form for the response function Riii(t) 

( ri(r) , r < 

^iii((r) ~ < ( 16 ) 

where the function ri(r) was computed by CD [S] while r(x) is given by 


'PedgeMdr . (17) 

Its small and large argument asymptotic can be easily obtained from the ones for Pedge(^") 


r[x) = e 

Jo 


in Eq. (10) with (3 = 1 and are given by 

1 


1 


r{x) ~ < 


4V^ a;3/2 


, X —)■ 0 , 


(18) 


ai 

— , X -)■ CX) , 

x^ 


where Oi is the amplitude in Eq. ( pTj ) for (3 = 1. The limit x —)■ 0 naturally matches the 
result of CD while the large-x behavior gives rise to a different algebraic decay. This late 
time regime III corresponds to the final stage of relaxation within the multidimensional 
energy landscape dominated by both the global minimum and the small-index saddle- 
points. 

Finally, for t and t' ~ the response and the correlation functions 

are characterized by scaling functions of the two variables t/N‘^3^ and t! whose 
asymptotic behaviors match smoothly with the various regimes described above. For 
instance, the response function i?(t, t'( takes the scaling form 


1 




.? = 


t' 


iV2/3 


,t>t\ 


(19) 


which is in general not stationary. Such a non-stationarity is usually interpreted as 
a form of ‘aging’ typical for complex systems [8]. However, in the limit tR' S> 1 


stationarity is restored so that Rftft') ~ r(t — t'), where r(x) is given in Eq. (17) 


3. Zero temperature dynamics 

The general solution of (|^, starting from a given initial condition sa( 0) at t = 0 for 
A G Sp(J) was found by CD and reads [8]: 


s\{t) = SA(0)exp(At)exp 


— / z{T)dT 


( 20 ) 


dt” exp [A(t — t”)] exp 


z{T')dT' 


>t" 


(hxin+un), 
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which will be the starting point of our analytical computations. In the following, we 
will focus on the uniform initial condition sa( 0) = 1, see Eq. 0- 


3.1. Lagrange multiplier z{t) and energy density e(t) 


At T = 0, the expression of z{t) is simply obtained from the spherical condition 
C{t,t) = 1. Using the definition of C{t,t) in Eq. @ together with the exact solution 


of the Langevin dynamics in Eq. (20), then setting hx{d) = ^x{d) = 0 for all 0 < f 
(and all A’s) one obtains: 


N = 


Vexp [-2 [ 
, L Jo 


{z{t) — X)dT 


After simple algebra the spherical condition (|21|) can be represented as: 

1 d 


z{t) = 


2 dt 




which by separating the contribution from the maximal eigenvalue An 
conveniently rewritten in the form 

1 d 


Z{t) = Amax + 2 ^ ’ 


( 21 ) 

( 22 ) 
can be 

(23) 


^2(A—Aniax)i 




From Eq. (23) it is obvious that 
lim z{t) = Amax • 


t^OO 


(24) 


The hnite time behavior of z{t) is therefore controlled by the random function gN{t) 
in (23). It is useful to write the latter function in terms of the density of near-extreme 
eigenvalues given in Eq. 0 as 

1 iV - 1 


9Nit) ~ Jf 


N 


-/lAf(t) , hN{t) = 


-2rt 


PDOs(r, N)dr . 


(25) 


From the results for Pdos(^) -^) [ 20 ] stated above in Eq. (10), we expect a very different 
behavior of hjq(t) for t ~ 0(1) (and more generally for times of order of f and 

for longer times t > 

Indeed, for t ~ 0{1), the integral in Eq. (25) is dominated by the region r ~ 0{1) 
where the DOS is (i) self-averaging and (ii) given by the shifted Wigner semi-circle law 
Pbuik(^") in Eq. (10). Therefore for this range of times we have: 

lim hjv(t) = hN(t) = h(t) = [ e“^''*pbuik(?") dr, for t ~ (9(1) . (26) 

V—^OO /n 


The integral over r can be performed explicitly [H] which yields [§] 


dr 




l‘{0 = / = ^A(4«) , 


( 27 ) 


§ Note that h{t) = e in the notation of Ref. [5]. 
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where Iy{x) is the modified Bessel function of index v. Eq. (27) implies that h{t) is a 
monotonically decreasing function of t whose asymptotic behaviors are given by 

1 — , t —y 0 , 


h{t) 


(28) 


4^^ 


, t > 1 . 


On the other hand, for t ~ one expects that the integral in Eq. (25) will be 

dominated by r ~ where the DOS behaves differently, see Eqs. (10) and (11). 

For such values of r one expects that the DOS is not a self-averaging quantity, which 
implies that hj\f{t), for t ~ is not self averaging either - and we have checked 

this statement numerically in Appendix A Finally, for very large times , t S> 
the function hN(t) is obviously dominated by the first gap between the two largest 
eigenvalues = Ai — A 2 , where we have used the ordering Amax = Ai > A 2 > • • ■ > Aat. 
We then arrive at the following asymptotic behaviors for hN{t): 

1 


/l7v(t) 


, 1 -C f < 77^/^ , 


4 ^^ 


—e f , 


(29) 


while in the crossover regime t ~ hj^^t) is a non-trivial random variable 

whose exact statistics is related to Pdos(^;^) at the edge and is presently not available 
analytically. Again we emphasize that hN{t) is self averaging only for t 

Before analyzing two-time quantities, it is already interesting to study the average 
value of the Lagrange multiplier z{t) which is related to the average energy density 
e{t) = i7[{si(f)}]. One has indeed the following relation (see also [8]): 


e(t) = 


z{t) 


A. 


h)v(t) 


41/{N-l) + hN{t) 


(30) 


For t < C>(iV2/3), one can replace hf^{t) by its large N value h{t) given in Eq. ( |27[ ), 
to obtain the result given in introduction in Eq. (13) with the function ei{t) simply 
given by 


1 hfit) 


, , 1 Wt) 

2t h(it) ’ 


(31) 


which has the asymptotic behavior given in Eq. (12). 


Let us now analyze this formula (30) in the limit t 3> where h]y{t) can be 


replaced by its asymptotic behavior given in (29). One finds: 

. 1 - 


Ae{t) ~ 


(32) 


The computation of this average requires the knowledge of the probability distribution 
function Pgap{g, N) of the first gap g, and in particular its small argument behavior since 
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the average in Eq. (32) will be dominated, for large t, by small gaps g. From Ref. 


one has (see also H) 

Pgv(s, ~ , Pg.p(x) « x** , for X -> 0 , (33) 

for general /9-ensemble. Hence, for fi = 1 which corresponds to onr case, one obtains 
that Ae(f) oc for t ^ as annonnced in Eq. (13). In the crossover 

regime t ~ Ae{t) is a non-trivial fnnction interpolating between the two 

limiting regimes (IT^, which is however hard to compute analytically. Thus we have 


demonstrated that the very late time dynamics is dominated by the hrst gap g between 
the two hrst eigenvalues of the coupling matrix. Interestingly, the same hrst gap also 
governs the huctuations of the overlap between two spin-conhgurations, see 


3.2. Response function 

In this section, we compute the response function R(t, t') dehned in Eq. 


From 


Eq. (20) we obtain an explicit expression for such a response function as (see Appendix 


[^for more details) 

R{tfi') = 



9N{t') 1 

9N{t) N 


E 


3 -(Amax-A)(t-t') 




l9N{t') 

g2Ai-(Amax-A)(t-t')- 

/ 9N{t) 



fi>t' . 


(34) 


(36) 


Starting from the above formula we will be able to demonstrate the existence of three 
diherent time regimes, as already mentioned in the introduction, which we analyze here 
in detail. 

• Regime I where 1 T <C and 1 f C in this case, 9 N{t') and 9 N{t) 


are self-averaging and one can easily check that the term in Eq. (35) is of order 0{1/N) 


while the term in Eq. (34) is actually of order 0{1). One obtains in this case |H]: 

3/4 


N- 




h ^ 


t' 


(36) 


where the function h{x) is given in Eq. (27). As already demonstrated in Ref. jH], the 
dynamics for t, t' -C still displays two diherent regimes of times, as can be seen 


from the structure of the response function in Eq. (36). The hrst regime corresponds to 
1 -C f, keeping f = r > 0 hxed. In such a regime the response function is 

time translationally invariant, Ri{t,t') ~ h(r/2), and so is the correlation function. In 
addition it was shown [8] that in that regime the response and the correlation function 
satisfy the huctuation-dissipation theorem (FDT). For that reason this regime is called 
a “quasi-equilibrium regime”. However, it is important to note that in such a regime I 
the system is not at equilibrium. Indeed, taking t, 3> 1 but keeping 0 < t'/t < 1 hxed 
one can see that the response function is clearly non-stationary. This is the so-called 
“aging regime” where, for t 3> t', one has: 

R{fit') -^-^ 




(37) 
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One can further show that in the aging regime the FDT does not hold [ 8 ] . 

• Regime II where 1 -C F <C and t 3> In this case hj^(t') is self-averaging 

and given by hN{t') l/{4:'/^)t' while hAr(f) 1/N to leading order. On the other 


hand, for t 3> one easily sees that the second term in (34) decreases exponentially 
with time t. Finally, as f — F 3> in this regime the sum over A 7 ^ A^ax in fhe hrst 


term in Eq. (34) is actually dominated by the hrst gap g at the edge of the spectrum 


of the coupling matrix J: 

Ritfi') 


^-9t 


After averaging over the distribution of the hrst gap g one obtains, for t ^ using 


Eq. (33): 


R{fit') 


Ai 


(39) 


This form (39) is however not very illuminating. Instead, it is more convenient to rewrite 


R{t,f) in the scaling form valid for 1 ^ F and t ~ 


R{t,t') 


-Jr 


t 


where fjx) has the following asymptotic behaviors 

, a; 1 , 

x~‘^ , X 3> 1 . 


(40) 


(41) 


Note that the small x behavior indeed matches the one given in (37), as it should, while 


the large x behavior is a new result. In this regime II only an aging regime is present 
as the time t is necessarily much larger than t'. 

• Regime III where both t' ^ A^^/s ^ ^ 77^/3^ ~ hM{t') ~ 


1/N while the term in (35) decays exponentially with t. Hence in this regime one hnds 
that the response function is stationary and given by 


R(f, t') ~ Riii(t 


F) = - ^ g—(-''max —A)(t—t') — 


PDOs(r)e-''(*-*')dr . (42) 


Hence one obtains the result announced in the introduction (15) 


3.3. Correlation function 

The correlation function C{t,t') can be computed straightforwardly at zero temperature 


T = 0 from Eqs. (20) and 

9 n 

C{t,f) = 


. The resulting expression is 
t + d' 


\JaN(t)afiit’ 


(43) 
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where the function gNit) is given in Eq. (25). As for the response function, there are 
three different regimes, in the {t,t') plane. 

• Regime I where 1 and 1 <C t In this case, gN{t), gN{t') and 

gN{(t 1')/2) are self-averaging and one finds 

, (44) 

which is the result of Ref. [H]. In this region of times the full function C{t,t') is self¬ 
averaging. 

• Regime II where 1 and t 3> In this case, only gN{t') is 


self-averaging. By analyzing the above formula (43), one obtains the scaling form 

t' 


Cifif) 


iV2/3 J 


fc 


t 

Ar2/3 


(45) 


The exact computation of the full scaling function fc{x) remains a challenge but its 
asymptotic behaviors are known: 

, X 1 , 

(46) 

1 , X > 1 . 

• Regime III where t' S> and t S> In this regime one has (y'Af(^) ~ 

gN(t') ~ gN{(t + t')/2) r'U 1/N. We then see that to the leading order the (zero- 
temperature) correlation function trivializes: 

C{fit') ~ 1 . (47) 

This very simple result reflects the fact that, in regime III, sx{fi) ~ s\{t) ~ 


4. Conclusion 

To conclude, we have studied the zero temperature relaxational dynamics of the spherical 
p = 2-spin glass model of size N when N is large but finite. We have identified a 
crossover time scale fcross ~ beyond which the relaxation occurs within the 

“well” close to the global minimum of the energy landscape and is dominated by saddle 
points with a few unstable directions. This last stage of the dynamics is controlled by 
the so called density of near-extreme eigenvalues pDos{r, N) ([^ of the coupling matrix, 
see e.g. Eq. (16), which was recently studied in the context of RMT HE EH]. 


The problem of the long time dynamics of the same model for finite N and finite 
temperature, T > 0, was addressed in Ref. [21]. The presence of finite temperature 
introduces a new relaxation time scale, which turns out to be the longest one in the 
system, related to the activation over the energy barrier separating the two ground 
states, that are related by spin inversion. The following simple argument proposed in 
m shows that the energy of the minimal height of this barrier is given simply by the 
difference A 2 — Ai. Indeed, denoting by |si) (resp., |s 2 )) the appropriately normalized 
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eigenvector of the interaction matrix corresponding to the largest eigenvalue Ai (resp., 
second largest A 2 ) we may consider the linear combination Is^) = cos0|si) + sin 0182 ). 
When the parameter 0 changes continuously between 0 = 0 and 0 = tt the vector |s6i) 
describes a path along our sphere connecting the two global minima in the associated 
energy landscape: |si) and —|si). As easily seen the energy along such a path is given 
by -E(0) = —y (Ai cos^ 6 ^ + A 2 sin^ 6 ') and is maximized at 6 ^ = 7r/2, that is exactly at 
the saddle-point with a single unstable direction and height given by A 2 . The barrier 
along this path is then obviously given by N(A 2 — Ai) and scales, in a large but finite 
system, like 

The authors of [23] further suggested that this property will be reflected in the 
asymptotic behavior of the finite temperature spin-spin autocorrelation function. We 
hope that the approach presented here can be useful to provide a detailed analysis of such 
an observable for finite N. We leave this challenging question for future investigations. 

In this paper, we restricted our study to the spherical p = 2-spin model. Of course, 
it is natural to wonder whether our results will, at least to some extent, hold for the 
spherical p-spin model with p > 2. This is a challenging open question as, for p > 2, the 
geometry of the energy landscape is much more complicated than for p = 2 , with many 
(exponentially with N) minima [251126] (for a rigorous proof see [27]). Similar questions 
can be naturally asked about relaxation in other mean-field models with finite number 
N of continuous degrees of freedom and exponentially many minima in the landscape, 
like those considered in [28l |29] . 
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Appendix A. The questions related to self-averaging 

In this appendix we study the property of being self-averaging for the density of near¬ 
extreme eigenvalues Pdos (’"5 N). For our purpose, the concrete question is whether gN{t), 
or equivalently hj^{t), is self-averaging. We have studied this question in detail in the 
case of GUE, as we have a full analytical description of the average DOS, Pdos(?", A7), 
in that case. The same qualitative behavior is supposed to be valid for GOE as well. 
The full computation of Pdos(^)-^) was performed for the Gaussian Unitary Ensemble 
in Ref. [19]. From it, one can infer the scaling form of Pdos(’" 5 A^) depending on the 
scaling of r with N. If r ~ (P(l), Pdos (^5 is given by a shifted Wigner semi-circle law, 
which was the case analyzed in [ 8 ]. On the other hand, if r ~ one expects 

that Pdos(^) will have a different form, which is in principle hard to compute. 
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Figure Al. Plot of Nhpi{t) as a function of for iV = 50 for GUE. In blue: 10 

different simulations for hpf{t). In red: /lAr(t) on 10^ simulations. The straight lines 
correspond to the algebraic behaviors and , with /3 = 2 here for GUE. 


We have tested numerically if the following statement holds: 


lim hAr(t) = / e '' pBosir, N)dr , 


N^OO 


(A.l) 


where means an average over the random matrix ensemble. Our analysis shows that 
h]s[{t) is indeed self-averaging for t ~ 0{1), for which the integral over r is dominated 
by the region r ~ 0{1) i.e. 


dr 




lim hi^{t) = h{t) = — e — r) = Ji(4t) , t ~ 0{1) (.A.2) 


N^OO 


2n 


2t 


As we show now, self-averaging, as stated in (A.l) does not hold when t rs-/ 
i.e. when the integral over r is dominated by r ~ Indeed, if self-averaging 


(A.l) held we would have 




,h{x)= e ^^^Pedgeiu)du . 


(A.3) 


In Fig. 


Al 


we show a plot of Nhj\f(t) as a function of for 10 different realizations 


of Af X Af GUE random matrices, with N = 50. We see clearly that the different curves 


do coincide for -Cl - corresponding to the bulk regime (A.3) but these curves 


clearly differ as soon as ~ 1, which is a clear indication that self-averaging does 

not hold in this regime. 


Appendix B. Computation of the response function R{t,t') 

In this appendix, we give the details of the computation leading to the expression for 


the response function R{t,t') given in Eq. (35). The starting point of this computation 
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is the expression for the response given in Eq. (|8|), for t > t'\ 

1 

N 




h=0 


where sx{t) evolves via the Langevin Eq. Q with ^\(t) = 0, for all A and t 
are interested here in the zero temperatnre T = 0 dynamics): 
dsx{t) 


dt 


= (A - z{t))sx{t) + hx{t) , 


(B.l) 


since we 


(B.2) 


starting from the initial condition sa(^ = 0) = 1 for all A. Here hx{t) represents an 
(inhnitesimal) external magnetic field and z{t) is a Lagrange mnltiplier imposing the 
spherical constraint J2x '^a(^) = Hence z{t) depends on hx(t') for all A and 0 < t' <t. 


The solntion of this Langevin eqnation is given in Eq. (20): 

/•i 


sx{t) = exp(Af) exp 


— / z{T)dT 


(B.3) 


+ / dt” exp [X{t — t”)] exp 


It" 


z{T')dT' 


hx{t'‘ 


From this expression we see that sx(t) depends both explicitly on hxifi) (see the second 
term) and implicitly throngh the Lagrange mnltiplier z{t). It is convenient to introduce 
the function 

ft 


9N,h{t) = exp 


dt'{z{t') - Amax) 


(B.4) 


such that gN,h=o{t) = gN{t) given in Eq. (23). In terms of gN,h{t) on has, from Eq. (B.3): 
6sx{t) 


dhxit') 


h=0 


gN{t) 


„(A-Amax)(t-t') _ 


2 [gN{t)f/^ 


S\-\^.Lt ^9N,h{t) 

Shx{L) 


h=0 


(B.5) 


h=0 


in Eq. (B.5) can be computed from the normalization condition 


The term 

Shx{t') 

X]a'^a(^) = Indeed, by expanding the latter relation, using Eq. (|B.3|), to lowest 
order in hx's one obtains 


E 


,2Ai-2 drzir) ^ drz{r) / 


z{T')dT' 


hx{t”) 


0{K 


= N 


(B.6) 


By differentiating this relation (B.6) with respect to hx{t') and setting hx{t”) = 0 for 
all 0 < f" < t, one finds [using the definition of gN,hit) in Eq. (B.4)], from the explicit 

= 0 . (B.7) 


solution in Eq. (B.3): 

1 dgN,h{t) 


[gN{t)Y 5hx{t') 

Therefore one has 

dgN,h{t) 


h=0 


g2(A Amax)i _j_ 2q- 


-(A,^ax-A)l-(A„,ax-A)p-t') V g^it') 

gN{t) 


dhxit') 


/i=0 
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By injecting this relation 
R(t,t') = 


into Eq. (B.5) one arrives at 


'gN(i') 1 


gN(t) N 


E' 


j-(Amax-A)(t-t') 


IgNit') 


gN{t) 




2\t 


. (B.9) 


which yields (by isolating the term A = Amax in the hrst sum) the expression given in 
the text in Eq. (35). 
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